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Abstract 

We report a study of magnetic dynamics in multiferroic hexagonal manganite HoMnOa by far- 
infrared spectroscopy. Low-temperature magnetic excitation spectrum of HoMnOs consists of 
magnetic-dipole transitions of Ho ions within the crystal-field split J = 8 manifold and of the 
triangular antiferromagnetic resonance of Mn ions. We determine the effective spin Hamiltonian 
for the Ho ion ground state. The magnetic-field splitting of the Mn antiferromagnetic resonance 
allows us to measure the magnetic exchange coupling between the rare-earth and Mn ions. 
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Superexchange is one of the main concepts of magnetism in solids. Recently, superex- 
change was found to mediate magnetoelectric (ME) coupling in multiferroics - materials with 
coexisting magnetic and ferroelectric ordersii^i^. ME coupling may find important technolog- 
ical applications, for example, in an electric- write magnetic-read memory element^. Multifer- 
roic hexagonal manganites display robust room-temperature ferroelectricity and remarkable 
ME behavior^ that allow technological exploitation^. In these multiferroics, the superex- 
change interaction between the rare-earth and Mn ions is responsible for the strong ME 
coupling that allows the control of ferromagnetism by an electric field^. The details of the 
rare-earth/Mn superexchange coupling have heretofore remained unknown. Here, we report 
the measurement of the rare-earth/Mn superexchange in multiferroic hexagonal HoMnOa 
via the detection of an antiferromagnetic resonance by far-infrared (far-IR) spectroscopy. 
These results demonstrate the ferromagnetic nature of the rare-earth/Mn exchange that 
enables the electric- field control of magnetism in HoMnOa. 

HoMnOs (HMO) crystallizes in a hexagonal lattice, space group PQ-^cm. The crystal 
structure consists of layers of corner- sharing trigonal MnOs bipyramids separated by lay- 
ers of Ho^"*" ions^. Ferroelectric polarization along the c axis {TpE = 875 K) results from 
electrostatic and size effects that lead to the buckling of MnOs bipyramids and the displace- 
ment of the Ho^"*" ions out of the a-b plane^'-Z. The magnetic structure is formed by Mn^+ 
ions in a two-dimensional triangular network in the a-b plane coupled by antiferromagnetic 
(AF) exchange and also by Ho^+ i^h) ions located at Ho(l) and Ho(2) lattice sites with 
point symmetries Csy and C3, respectively^"^. The AF exchange between Mn^"*" magnetic 
moments leads to their 120° ordering at T/v = 72 K— lii, and easy plane anisotropy confines 
the Mn spins to the a-b plane. Two in-plane spin reorientation transitions occur at temper- 
atures TsRi^38 K and Tsr2^5 K associated with changes of magnetic symmetry to PQ'-^cm' 
and to PG^cm, respectively^>i2i>i^. The spin reorientation at Tsr2^5 K is accompanied by 
antiferromagnetic ordering of the Ho^+ ions located at the Ho (2) lattice sites^. Several man- 
ifestations of ME coupling in HMO have been documented, one of the most prominent being 
the presence of the reentrant phase with a strong magnetodielectric response^»^>^»^>^. 

The HMO single crystals were grown using an optical floating zone furnace-^. The ma- 
terial's magnetic phase diagram^^ at temperatures above 5 K displays the reentrant phase 
discovered by Lorenz et alM. Polarized far-IR transmission of a 0.4 mm thick, (110) ori- 
ented crystal was measured as a function of temperature and applied magnetic field using a 
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FIG. 1: Power transmission spectra of HoMnOs with polarized far-IR light incident along the 
crystal's [110] direction at T=10 K. Arrows indicate magnetic absorption lines. The dotted line 
illustrates the simulation of the transmission spectrum using Eqs. ([T]) and ([2]) and the parameters 
eo = 9.07, a = 0.9 cm^^, with magnetic resonances centered at 27.5, 43.5, and 52.1 cm~^. 

Bruker 66 Fourier-transform spectrometer coupled to an 8 Tesla split-coil superconducting 
magnetji^. Static magnetic field B was applied both parallel and perpendicular to the c axis. 

Figure [1] shows far-IR power transmission spectra of the HMO crystal at low temperatures 
and zero magnetic field for two polarizations of the light incident along the [110] direction - 
the magnetic /i-field of the lightwave parallel {h\\c) and perpendicular (/i±c) to the c axis. 
The spectra in Fig. [1] are normalized to the transmission through an empty aperture of the 
same size as the sample. The fringes with the period of ~ 3cm~^ seen in the spectra are due to 
multiple reflections of light within the sample. The transmission minima indicated by arrows 
correspond to magnetic resonance frequencies. The broad absorption in both polarizations 
at frequencies above 70 cm~^ is of non- magnetic origin, as it exhibits no magnetic field 
dependence. The dotted line in Fig. [1] shows the simulated transmission, assuming that the 
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magnetic and dielectric response of HMO consists of a collection of lorentzian oscillators: 



where uij, 7^, and fj are frequency, relaxation rate, and oscillator strength of magnetic 
resonances, eo and a are constants, and c^o, 7o, and fo are introduced to account for the 
non-magnetic absorption above 70 cm~^. 

Color-coded transmission maps in Fig. [2] illustrate the magnetic field dependence of the 
resonance frequencies at T=10 K. To create the maps, we recorded the transmission spectra 
in 1 T steps between and 8 T. The darker color in Fig. [2] indicates lower transmission 
and the white color indicates where the transmission is the highest; the symbols specify 
the positions of magnetic resonances extracted from the numerical simulation of the far-IR 
transmission. In the remainder of this Letter, we show that the magnetic field dependence of 
the resonance frequencies carries an unambiguous signature of the superexchange interaction 
between magnetic Ho'^"^ and Mn^^ ions that allows us to quantify the interaction. 

We attribute the absorption lines labeled 1-4 in Figs. [1] and [2] to crystal field excitations 
of Ho'^^ ions. The crystal fields of C^^ and C3 symmetry have qualitatively similar effects 
on the Ho ion. According to the work of Elliott and Stevensi^, the crystal field potential 
VcF can be written in terms of spin operators J, J^, and J-t for the manifold of states of 
constant J. The diagonal part of Vcf depends only on J and Jz- The off-diagonal part of 
VcF has non-zero matrix elements only between the states m and m' that differ by Am = 6, 
m' = m ± 6, in the field of C3 symmetry and by Am = 3, 6, m' = m ± 3, 6, in the field 
of C^y symmetry. In zero applied magnetic field, the degenerate J = 8 manifold splits 
into a set of singlet and doublet states of the form |±6, ±3, 0) and |±8, ±5, ±2, =i=l, =i=4, =i=7), 
respectively. We restrict our discussion to a qualitative description of the observed crystal 
field transitions, as the frequency range of 20 — 80cm~^ in our far-IR measurements does not 
allow us to characterize all magnetic-dipole active crystal field transitions and determine 
the crystal field parameters. For example, the 12.5 cm~^ transition observed by inelastic 
neutron scatteringi& is outside of our accessible frequency range. 

The magnetic-dipole selection rule for the (/i_Lc) polarization allows only transitions with 
Am = ±1. This rule forbids transitions between singlets, but allows transitions between 



e(u;) = eQ + i4:Ha/u + 



fo 
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FIG. 2: Transmission maps for the (/i_Lc) polarization where darker color indicates lower transmis- 
sion. Left panel - static field B\\c, right panel - i?_Lc. The bars above the panels set the absolute 
transmission scale. Symbols represent magnetic resonance positions extracted from fitting the 
measured transmission using Eqs. ([T]) and Resonances are labeled by the same numbers used 
in Fig. [TJ The inset of the left panel is a scheme of the transition 2, which happens from the 
ground state (GS) doublet. The final state of the transition is assumed to be a singlet. In low 
magnetic field, the GS doublet splits and two resonance lines are observed as illustrated by the two 
vertical arrows on the left. In high magnetic field, the spitting of the GS becomes large enough 
for a thermal depopulation of the upper branch. Only one resonance line is observed, as the single 
vertical arrow on the right illustrates. 

singlets and doublets and between different doublets. Transitions 1 and 2 are strong in the 
(/i-Lc) polarization (Fig. [2]) but are very faint in the {h\\c) polarization, which suggests that 
they are allowed in the former, but forbidden in the latter, polarization. Since we observe 
transition 2 in the highest applied magnetic fields, it must be the transition from the ground 
state. The final state of the transition could be either a doublet or a singlet. With a magnetic 
field B applied along the c axis (left panel of Fig. [2] and its inset) the transition splits linearly, 
a characteristic of the Zeeman splitting of a doublet. The lower branch of the Zeeman-split 
transition disappears when the applied field increases beyond 3 T, which is indicative of a 
thermal depopulation of the upper branch of the Zeeman-split ground state doublet. The 
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slope of the upper branch of transition 2 allows an estimation of the ^f-factor in the effective 
spin Hamiltonian Tie// = gf^'sBS, S = 1/2, that describes the ground state. Assuming that 
the final state of the transition is a singlet with energy independent of the applied field (at 
sufficiently low fields), we get ~ 8.8. The magnetic field B perpendicular to the c axis 
produces no Zeeman splitting of the ground state doublet. Therefore, neglecting the mixing 
of the ground state with higher energy levels at sufficiently low field B, we set Qx = Qy = 0. 
The upper branch of transition 2 is further split by about 5 cm~^ in magnetic fields along 
the c axis, which we attribute to hyperfine interactions. 




FIG. 3: Field dependence of triangular AFMR frequencies at T=10 K. The measured splitting of 
the resonance by magnetic field is about twice the splitting calculated using the free energy of Eq. 

We now turn to the magnetic response of Mn ions in HMO and assign the sharp absorption 
line at 43.4 cm^^ in the (/i_Lc) polarization to the antiferromagnetic resonance (AFMR) of 
the triangular lattice of Mn magnetic moments (transition 5 in Figs. [1], [2]). This assignment is 
supported by the observation of the triangular AFMR at the same frequency (~ 43cm~^) and 
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light polarization in hexagonal YMnOjr^, a related compound with the same Mn magnetic 
structure and similar exchange and anisotropy parameters'^. The AFMR corresponds to the 
gapped k = magnon detected by inelastic neutron scattering at 47 cm~^ in HMO^^. The 
similarity between the k = magnon frequency and the position of the AFMR in our data 
further supports the assignment of the resonance. The magnetic field dependence of the 
AFMR agrees qualitatively with the expected linear splitting induced by the field applied 
along the c axis (-B||c) (left panel of Fig. [2]). A slight increase in the AFMR frequency and 
no splitting of the resonance are observed in fields up to 5 T applied perpendicular to the 
c axis (-B_Lc). In fields of 5 T and higher, the AFMR line is hard to identify, as it becomes 
a part of the broad absorption band due to the crystal field levels of Ho ions (right panel 
of Fig. [2]). The same qualitative behaviour, i.e., linear splitting with B\\c and no splitting 
with i?±c, in fields up to 5.4 T was recorded for the AFMR in YMnO^^. 

To quantify the magnetic-field behaviour of the AFMR, we use the following free energy: 

F = Xj^Mi - Mj + Kj^iM^y - bJ2m:, (3) 

ij^j i i 

where Mi, i = 1, 2, 3, are the three sublattice magnetizations within the a — b plane, A 
describes the antiferromagnetic exchange between the sublattices, and K represents the 
easy-plane anisotropy. The B = excitation spectrum of the free energy ([3]) consists of a 
zero-frequency mode corresponding to global rotations of the sublattices about the c axis 
and a doubly-degenerate gapped mode whose degeneracy is lifted by the application of i?||c. 
The analytical expression for the frequency of the gapped mode in B\\c was given by Palme 
et al^: 

2 _ «^ b{a-b) 2 
bB 



'B^b{b-2a) + 2ab{a + bY, (4) 



2(a + 6)2 

where a = 2Hd, b = 3Hex, the exchange field Hex = AMq, the anisotropy field = KMq, 
and Mq is the sublattice magnetization. Equation ([3]) is a molecular-field version of the 
Hamiltonian used by Vajk et al to describe the spinwave dispersion in HMO^^. In this 
picture. Hex = 35*7 = 127 T and Hd = SD = 5.75 T, where J and D are the exchange and 
anisotropy parameters measured by inelastic neutron scattering and 5* is the total spin of 
Mn ions. The value of H^, was lowered slightly to reproduce the zero-field AFMR frequency 
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observed in our far-IR study. The field dependence of tlie AFMR calculated using Eq. 
(jl]) and the (/-factor g = 2 is shown in Fig. [31 The calculated frequencies largely disagree 
with the measured ones as the measured splitting is about twice the calculated splitting. For 
comparison, the splitting of the AFMR in YMnOs is in good agreement with our calculation: 
the measured splitting of 4.8 cm^^ at ~ 5.4 T {g = 1.9 ± 0.1)^^ agrees with the calculated 
value of 4.9 cm~^, as determined from exchange and anisotropy fields measured by inelastic 
neutron scattering2S in YMnOg^. 

What causes the large discrepancy between the measured and calculated splitting of 
the AFMR line in HMO? Our calculation describes well the AFMR in YMnOs, a material 
that is largely similar to HMO, but with one important difference: the Y'^'^ ions are not 
magnetic. Therefore, we conclude that the magnetic Ho ions in HMO are responsible for the 
discrepancy. The Mn ions in HMO interact with the surrounding Ho ions via an exchange 
coupling of the form J-^jf^" ■ Sj^"^, where S^'^ and J^" are the angular momenta of Mn 
and Ho ions. Using the equivalence of matrix elements gjJz = gzSz, where gj = 5/4 is the 
Lande factor, gz and 5*2 are the effective g-factor and spin for the Ho ion ground state, we 
can write the exchange Hamiltonian as^® Hex = JxSxS^^^ + JySyS^'^ + J'zSzS^'^"', where 
J^z = {gz/gj)J'ij, etc. We found earlier that gz ~ 8.8 and g^ = gy = 0, which reduces 
the Hamiltonian to Tiex = J'zSzS^^"'. In magnetic field B applied along the c axis, we can 
replace Sz by its thermal average (Sz)= xB / gzfJ-B and write Hex = i>JzX)/i9zfJ'B)BSf^'^, 



where x is the magnetic susceptibility of Ho ions. The last expression shows that the Ho-Mn 
(HM) exchange interaction is equivalent to an effective magnetic field acting on Mn spins. 
To describe the effect of the HM exchange on the AFMR, we add the term 



where N=Q is the number of Ho ions neighboring a Mn ion, to the free energy of the Eq. 
([3]). The total magnetic field acting on Mn ions is then the sum of the external static field 
and the effective field of Ho ions, = (1 — Xhm)B, and can be larger or smaller than the 
external field, thus enhancing or mitigating its effect on the AFMR. Comparing the measured 
AFMR splitting and the splitting calculated using Eq. Q, we find that Xhai = —1, which 
corresponds to ferromagnetic HM coupling. 

We can account for the two kinds of Ho lattice sites (Ho(l) and Ho(2)) by writing J'z = 
^Jl + ^Jz')l^ 1 where Jl and are the HM exchange constants for each site. Each Mn 





(5) 
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ion has two Ho(l) ions and four Ho(2) ions as Ho nearest neighbors. The two Ho lattice sites 
are distinct because of the ferroelectric displacements. Such displacements are significantly 
smaller that interionic distances in HoMnOs, and the approximation Jl ^ ^ is 
expected to be a good one. Otherwise, we must use the above combination of Jl and J'^ 
in the definition of \hm in Eq- dS])- 

Before we conclude, we consider the possible effect of phonon coupling on the frequencies 
of the studied magnetic excitations. The most prominent signature of such coupling is the 
anti-crossing of the coupled magnetic and lattice modes. Such anti-crossing is expected to 
be strongest when the frequencies of the two modes are similar, as shown on the right panel 
of Fig. [21 where the strong magnetic absorption line overlaps at high magnetic field with the 
non- magnetic absorption at 70cm~^. At all magnetic fields, we see no measurable changes to 
the 70 cm~^ excitation. Thus, we conclude that we detect no measurable coupling between 
magnetic and lattice vibrations. 

The measured ferromagnetic HM exchange plays a significant role in all aspects of the 
physics of HoMnOa. It causes the electric-field induced ferromagnetic ordering of Ho mag- 
netic moments along the c axis^ and may also be responsible for the multitude of low- 
temperature (T<8 K) magnetic phases and the two Mn spin reorientation transitions in 
HMO^°. The reentrant phase that is induced in HMO by magnetic field displays a strong 
magnetodielectric effect and was explained by the formation of magnetic domain walls dur- 
ing spin reorientation, which reduces the local magnetic symmetry and allows the coupling 
between magnetic and ferroelectric polarizationsi^*^. Theoretical descriptions of such mag- 
netoelectricity must include domain wall contributions to the free energy from lattice distor- 
tions and magnetic superexchange interactions. Such descriptions are incomplete without 
the consideration of the magnetic exchange between the Ho and Mn ions. Our measurement 
of the ferromagnetic HM exchange provides an important contribution to the modeling of 
the interplay between magnetism and ferroelectricity in HoMnOs. 

We thank Darryl Smith for useful theoretical discussions and Jonathan Gigax for the 
simulation of the Ho crystal field eigenstates. The work at LANL was supported by the 
LDRD program and the Center for Integrated Nanotechnologies. The work at UCSD was 



9 



supported by NSF DMR 0705171. 



Electronic address: diyar@lanl.gov 



1 N. Hur, S. Park, P. A. Sharma, J. S. Ahn, S. Guha, and S.-W. Cheong, Nature 429, 392 (2004). 
^ T. Lottermoser, T. Lonkai, U. Amann, D. Hohlwein, J. Ihringer, and M. Fiebig, Nature 430, 
541 (2004). 

^ S.-W. Cheong and M. Mostovoy Nature Mat. 6, 13 (2007). 

4 J. Scott, Science 315, 954 (2007). 

5 M. Fiebig, T. Lottermoser, D. Frohlich, A. Goltsev, and R. Pisarev, Nature 419, 818 (2002). 
^ N. Fujimura, T. Ishida, T. Yoshimura, and T. Ito, Appl. Phys. Lett. 69, 1011 (1996). 

B. B. V. Aken, T. T. M. Palstra, A. Filippetti, and N. A. Spaldin, Nature Mat. 3, 164 (2004). 

8 A. Munoz, J. A. Alonso, M. J. Martinez-Lope, M. T. Casais, J. L. Martinez, and M. T. 
Fernandez-Diaz, Chem. Mater. 13, 1497 (2001). 

9 M. N. Iliev, H.-G. Lee, V. N. Popov, M. V. Abrashev, A. Hamed, R. L. Meng, and C. W. Chu, 
Phys. Rev. B 56, 2488 (1997). 

^° O. P. Vajk, M. Kenzelmann, J. W. Lynn, S. B. Kim, and S.-W. Cheong, Phys. Rev. Lett. 94, 
087601 (2005). 

11 N. Hur, I. K. Jeong, M. F. Hundley, S. B. Kim, and S.-W. Cheong, cond-mat:0805.3289 (2008). 

M. Fiebig, D. Frohlich, K. Kohn, S. Leute, T. Lottermoser, V. V. Pavlov, and R. V. Pisarev, 

Phys. Rev. Lett. 84, 5620 (2000). 
13 T. Lottermoser and M. Fiebig, Phys. Rev. B 70, 220407(R) (2004). 

1^ B. Lorenz, A. P. Litvinchuk, M. M. Gospodinov, and C. W. Chu, Phys. Rev. Lett. 92, 087204 
(2004). 

1^ C. dela Cruz, F. Yen, B. Lorenz, Y. Q. Wang, Y. Y. Sun, M. M. Gospodinov, and C. W. Chu, 

Phys. Rev. B 71, 060407(R) (2005). 
16 F. Yen, C. R. dela Cruz, B. Lorenz, Y. Y. Sun, Y. Q. Wang, M. M. Gospodinov, and C. W. 

Chu, Phys. Rev. B 71, 180407(R) (2005). 
i'^ W. J. Padilla, Z. Q. Li, K. S. Burch, Y. S. Lee, and D. N. Basov, Rev. Sci. Instr. 75, 4710 

(2004). 

1^ A. Abragam and B. Bleaney, Electron paramagnetic resonance of transition ions (Clarendon 

10 



Press, Oxford, 1970). 

T. Penney, P. Berger, and K. Kritiyakirana, J. Appl. Phys. 40, 1234 (1969). 

T. J. Sato, S.-H. Lee, T. Katsufuji, M. Masaki, S. Park, J. R. D. Copley, and H. Takagi, Phys. 

Rev. B 68, 014432 (2003). 

W. Palme, F. Mertens, O. Born, B. Luthi, and U. Schotte, Solid State Comm. 76, 873 (1990). 
Sato et ali^ find that the chemical triangular lattice of Mn ions is slightly trimerized with 
exchange strength Ji between the nearest neighbors belonging to the same trimer and exchange 
strength J2 between the ions belonging to neighboring trimers. The zero wave-vector magnon 
spectrum of the trimerized chemical lattice is calculated using Eq. ^ and the substitution 
(Ji + 2J2)/3. 



11 



